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We present a simple formalism for the evaluation of the Casimir energy for two spheres and a 
sphere and a plane, in case of a scalar fluctuating field, valid at any separations. We compare the 
exact results with various approximation schemes and establish when such schemes become useful. 
The formalism can be easily extended to any number of spheres and/or planes in three or arbitrary 
dimensions, with a variety of boundary conditions or nonoverlapping potentials/nonideal reflectors. 
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I. INTRODUCTION 

In 1948 the Dutch physicist Hcndrik Casimir predicted 
the existence of a very peculiar effect, the attraction 
between two metallic uncharged parallel plates in vac- 
uum Jj. The existence of such an attraction has been 
confirmed experimentally with high accuracy only re- 
cently 0, 01 ■ However, nearly all modern experiments 
(the noted exceptions are the two-cylinder work of Ref. 4] 
and the two-plate experiment of Ref. 0) study the at- 
traction between a metallic sphere and a metallic plate 
which are much simpler to align than two plates, but 
much harder to calculate. In fact, with the exception of 
the proximity- force approximation ^ i2j i which is only 
applicable for vanishing separation, there does not ex- 
ist a theoretical prediction for the Casimir energy of the 
sphere-plate system as function of the distance. 

The origin of this attractive force can be traced back 
to the modification in the spectrum of zero point fluctua- 
tions of the electromagnetic field when the separated mir- 
rors are brought into close distance. Similar phenomena 
are expected to exist for various other (typically bosonic) 
fields 13 and the corresponding forces are referred to 
as Casimir or fluctuation-induced interactions. A related 
interaction arises when the space is flUed with fermions, 
which is particularly relevant to the physics of neutron 
stars US El in mill Hi and quark gluon plasma [l^. 

Since the Casimir effect between a sphere and a plate 
is the experimentally interesting but theoretically dif- 
ficult case, because of the electromagnetic nature of 
the fluctuating fields, the corresponding Casimir ef- 
fect for a real scalar field between two spheres or one 
sphere and a plate came into the focus of theoretical re- 
search [HllllllllSElElEllllEl- Therefore we 
will here focus on the exact and semiclassical calculation 
of this scalar Casimir effect between spheres or spheres 
and plates in three spatial dimensions, where we assume 
Dirichlet boundary conditions on the spheres and plates. 
This scenario can be trivially extended to the case of Neu- 
mann or other boundary conditions, to two-dimensional 
systems of disks and/or lines and to analogous systems 
in arbitrary dimensions. We will show that the Casimir 



force calculation for the two-sphere case incorporates the 
sphere-plate geometry as special case. Note that, with 
the exception of the numerical results of Refs. [l^jllol (see 
also |23)> which do not yet extend to large separations, 
no exact result exists in the literature for the Casimir ef- 
fect for a real scalar field under Dirichlet boundary con- 
ditions, neither between a sphere and a plate nor be- 
tween two spheres. For small separations between the 
sphere and the plate, the above-mentioned proximity- 
force approximation [y, Q can be applied. Its justifica- 
tion at small separations has been provided by many au- 
thors using various theoretical techniques. For instance, 
in Refs. 0, 0| semiclassical methods in the framework 
of the Gutzwiller trace formula |26| | have been used, in 
Ref. '27*1 the proximity-force approximation for the elec- 
tromagnetic case has been derived from the multiple scat- 
tering expansion of Ref. p3 |. in Refs. |l9l |20| the world- 
line approach in the framework of the Feynman path in- 
tegral has been applied, and in Refs. 22, 23, 24, 25] a ray- 
dynamical approach in terms of optical paths (i.e. closed 
but not necessarily periodic orbits) has been employed. 

Here, we will present an evaluation of the scalar 
Casimir problem that is based on quantum mechanics, 
without any semiclassical approximation made before- 
hand. It utilizes the Krein formula [29l l30l | as a bridge 
between the spectral density on the one hand and the 
problem of scattering of a point particle between spheres 
in three dimensions (or disks in two dimensions) on the 
other hand. It has to be emphasized that in this case 
the Casimir calculation is not plagued by the removal of 
diverging ultra-violet contributions. This is related to 
the fact that the Krein formula is exact, and that the de- 
terminants of the S-matrix and the corresponding inverse 
multi-scattering matrix for a system of N nonoverlapping 
spheres (or disks) , which both are manifestly known, are 
finite [U m Ha. In this work we do not consider the 
material dependent stress of, e.g., the deformation of a 
single spherical shell which was discussed in Ref. • 

The paper is organized as follows: first, we formulate 
the problem in terms of the density of states of the scalar 
field and relate it to the S-matrix of the system using the 
Krein formula. In Section IIIII we focus on the particu- 
lar realization of the Casimir effect for the sphere-sphere 
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and the sphere-plate system. Section IIVI is devoted to 
investigations of the large-distance limit between spheres 
where the asymptotic expressions are derived and dis- 
cussed. Section |3 discusses the link between the pre- 
sented approach and the semiclassical methods based on 
the Gutzwiller trace formula. Finally, the numerical re- 
sults, approximate expressions and conclusions are pre- 
sented in Sections I VII and rVIII For the sake of complete- 
ness, the derivation of proximity-force approximation for 
the sphere-sphere and sphere-plate system is given in Ap- 
pendix ^ and a comparison to the two-dimensional two- 
disk and disk- line systems can be found in Appendix IbI 



II. THE MODIFIED KREIN FORMULA 

Our main goal is to reach a qualitative understand- 
ing of the scalar Casimir energy at zero temperature 
in the case of more complicated geometries than the 
original two-plate system. For that purpose let us con- 
sider the fluctuating real scalar field between N nonover- 
lapping, nontouching, impenetrable spheres of radii 
{i = 1, . . . , N). It is assumed that that the scalar field is 
noninteracting and is subject to Dirichlet boundary con- 
ditions on the surfaces of the spheres. The spheres are 
positioned at fixed relative distances r.y — Lij + + aj > 
tti + ttj between their centers; is then the shortest rela- 
tive distance between their surfaces, and fij is the center- 
to-center distance vector, which includes also the infor- 
mation about the spatial orientation. In order to calcu- 
late the Casimir energy we shall represent the smoothed 
bosonic density of states of the scalar field as a function 
of the energy e (smoothing is over an energy interval Ae 
larger than the level spacing in the big volume V of the 
entire system): 

N 

g{e, {ai}, {nj}) = go{e) + ^ gw{e, fli) 

+ gcie,{ai},{nj}) , (1) 

where ^(e, {a;}, {r^ }) is the total density of states of 
the scalar field, gois) is the density of states in the ab- 
sence of all scatterers, and gwi^^ on) is the correction to 
the density of states arising from the presence of one 
sphere (sphere i). Clearly 'Yl!i=i 9w ^ cn) is the correc- 
tion due to the N spheres infinitely far apart from each 
other that sums up the excluded volume effects, surface 
contributions and Friedel oscillations caused by each of 
the obstacles separately. Finally gc{£, {oj}; {^y})) is the 
remaining part, which is of central interest to us here. It 
vanishes in the limit of infinitely separated scatterers and 
is the only term in the density of states which reflects the 
relative geometry-dependence of the problem. Only this 
term contributes to the Casimir energy. 

Strictly speaking, only the smoothed level densities 
gwi.^) and gci^) are flnite, whereas the level densities 
g{e) and go{e) are infinite, as they are proportional to 



the volume V of the entire space. This redundant diver- 
gence can be handled easily by considering the smoothed 
quantities first in a very big box, the volume of which is 
subsequently taken to infinity 31]. 

Now we will use the Krein formula [s^l which 
provides a link between the (A^-body) scattering matrix 
Sn {s) of a point-particle scattering off A^ spheres and the 
change in the density of states due to the presence of A^ 
scatterers, namely 



5g{e, {ai}, {f^ }) = g{£, {aJ, {nj}) - 50(e) 
= 2^ ^lndetS'7v(£, {aj,{fy}) . 



(2) 



Note that In det Sn (e, {oi}, {^ij}) /2i is nothing else than 
the total phase shift of the scattering problem. The 
geometry-dependent part of the density of states can now 
be extracted from the genuine multi-scattering determi- 
nant. In this way the calculation is mapped onto a quan- 
tum mechanical billiard problem that classically corre- 
sponds to the hyperbolic (or even chaotic) point-particle 
scattering off A^ spheres ISal (o r A^ disks in two dimen- 
sions llillillllsilllliflllillil). 

As shown in Refs. |^ Is^ . l33| , the determinant of the 
A^-scatterer S-matrix, S'at (e, {a^}, {ry}) , factorizes into 
the product of the determinants of the single-scatterer 
S-matrices and the ratio of the determinants of the in- 
verse multi-scattering matrix 13 ll 133, l3J of Korringa- 
Kohn-Rostoker (KKR) type gl H El El M{k) = 
M (k , {ai} , {rij}) in the complex wave number (k = \k\) 
plane: 



N 



detS'Ar(£, {ai},{fy }) = < ]^detS'i(e, Oi) 



detM(r) 
detM(fc) 



(3) 

The formula Q holds in the case when the scattering 
modes are free massless fields as well as in the case of free 
nonrelativistic fields with a mass m. Both cases imply 
different energy dispersion relations, e = fku = hck in 
the massless case or e = h^k'^ / {2m) in the nonrelativistic 
scenario, respectively. 

Although the involved matrices are infinite- 
dimensional, all determinants are well-defined, as 
long as the number of spheres is finite and the spheres 
do neither overlap nor touch. This follows from the 
trace-class property [i^, E3 of the matrices Sn — 1, 
Si-t and M - 1 which was shown in Refs. |3ll32[l33l|. ^ 

Inserting the exact expression Q into the original 
Krein formula ((SJ, using the decomposition |^ and iden- 
tifying the "Weyl-type" density of states with the phase 



^ Trace-class operators (or matrices) are those, in general, non- 
Hermitian operators (matrices) of a separable Hilbert space 
which have an absolutely convergent trace in every orthonor- 
mal basis. Especially the determinant det(]l + zA) exists and is 
an entire function of 2, if is trace-class. 
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shift of the corresponding single scatterer 

gwi£,ai) = lndetS'i(£,ai) . 

2ni de 



(4) 



one finds a new Krein-type exact formula which di- 
rectly links the geometry-dependent part of the density 
of states with the inverse multi-scattering matrix 

9c{£,W^}AnJ}) = ^^ImlndetM(fc(£),{aJ,{r;.,}) 



^fc{£,{a^},{nj}) = --ImlndetM(fc(e),{aO,{r=;j}) 

(5) 

for the integrated geometry-dependent part of the density 
of states 

^^c{e,{a^},{nJ})= f de'gc{e',{a,},{r,,}). (6) 
Jo 

The exact formula © is the central expression of this 
paper. The Casimir energy itself follows via the integral 



£c = 



/ de^egcie) = -5 / deJVc{e) 
Jo Jo 

delm lndetM(fc(e)) 



1 

2^ JO 

1 
47ri 



oo(l+iO+) 



de Indct Af(fc(£)) 

D(l-iO+) ■ 

de lndetA/(fc*(£))'f 



(7) 



In the massless case e — Hck, this expression can be Wick- 
rotated (i.e., k ^ iki for the first term and fc — > —ik^ 
for the second term of the last relation) to give ^ 



he 



dfc4 Indet Miik^) 



(8) 



since det M{k) = detM((-fc)*)t and therefore 
detM(iA:4) = detM(iA:4)1' if k4 real. ^ 



Using the exp licit formulas for the KKR-type matrix 
from Ref. |3ll 133 . |33| . one can compute numerically 
£c for various arrangements of hard spherical (or circu- 
lar in 2D) scatterers. For a point-particle scattering off 
N nonoverlapping nontouching spheres (under Dirichlet 
boundary conditions) 



ME 



Irn.l' m' 



<5ii'^w + (1 - <5^^ )V47rf 



;2m-|-/ 



V(2; + 1)(2Z' + 1) (^) 



]i{kaj) 



X \/2/" + li'" 



I" 



V I 



0/ \m-m" m" 



(9) 



is the inverse multi-scattering matrix 33]. Here — 
1,2,...,N are the labels of the N spheres, I, I', I" = 
0,1,2,... are the angular momentum quantum numbers, 
and m,m',m" the pertinent magnetic quantum num- 
bers. DJj^, ^„(j, j') is a Wigner rotation matrix which 
transforms the local coordinate system from sphere j' to 
the one of sphere j, j/(fcr) and ]if\kr) are the spher- 
ical Bessel and Hankel functions of first kind respec- 
tively, YJI)"^ i^jf) ^ spherical harmonic (where fj^^ 
is the unit vector, measured in the local frame of sphere 
j, pointing from sphere j to sphere j'), and the 3j- 
symbols 49] result from the angular momentum cou- 
pling. 

By definition, the inverse multi-scattering matrix in- 
corporates the pruning rule that two successive scatter- 
ings have to take place at different scatterers (see the 
(1 — 5^^ ) term in Eq. This alternating pattern 

between a single-scatterer T-matrix and the successive 
propagation to a new scatterer IS lisf distinguishes 
the KKR-type method from the multiple scattering ex- 
pansion of Refs. [2112^. 



^ The Wick-rotations are allowed, since detM(fc) has poles in the 
lower complex fc- plane on l y, w hereas det M(fc*)t has poles in the 
upper half-plane l3lt l32l l33l . Furthermore, the integrals over 
the circular arcs vanish, since Indet M(fe) and Indet M(fc*)t are 
exponentially suppressed for Imk — > ±00, respectively; see, e.g., 
the semiclassical expression l4Ui . 

^ For the same reason, one can show the corollary that all the 
corresponding integrals over odd powers of k have to vanish: 

he 

= — / dfcfc2"+ilm Indet M(fc) 
27r Jo 

= i(-l)" — / dfc4 fc?"+M Indet M(ife4) - Indet M(iA:4)'l' . 

47r Jo L J 

Thus, the Casimir energy over modes with a nonrelativistic dis- 
persion e = ft^fc^/(2m) integrates to zero, unless there exists a 
finite upper integration limit, as e.g. the Fermi m omentum in 
the fermionic Casimir effect studied in Ref. Il3tll5l . 



III. THE TWO-SPHERE AND THE 
SPHERE-PLATE PROBLEM 



In the case of two spheres the system possesses a con- 
tinuous symmetry, i.e. C'oov in crystallography group 
theory notation ,5(|] , associated with rotations with re- 
spect to the axis joining the centers of the spheres (and an 
additional reflection symmetry with respect to any plane 
containing this symmetry axis). As a consequence the 
KKR-matrix is separable with respect to the magnetic 
quantum number m (in fact, it depends on its modulus 
|m| only) |33|. In the global domain, it is given as 



Im.l'' 



21 (m) 
W 



Sw 



(10) 
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with 



AJJ 



(-1)' i' 



Y^f{2l" + l) 



I" I V 



I" I V 




V(2; + 1)(2Z' + 1) 
4V(fcr). 



m —m 



(11) 



If the spheres have moreover the same radius ai = 02 = 
a, there exist also a two- fold reflection symmetry with 
respect to the vertical symmetry plane. This additional 
symmetry makes the total symmetry of the system to 
be Dooh in the crystallography group notation 50] 
which is a simply product of the Coov group and the 
inversion (and rotation by tt) with respect to the point 
of intersection between the symmetry axis and vertical 
symmetry plane. 

Therefore the global domain of the two-sphere system 
can be split into two half-domains, separated by this 
plane, see Fig.^ and all the (scattering) wave functions 
can decomposed into symmetric and antisymmetric ones 
with respect to the vertical symmetry plane. The sym- 
metric wave functions are subject to Neumann boundary 
conditions and the antisymmetric are subject to Dirichlet 
boundary conditions on this symmetry plane. Thus there 
exist two KKR matrices in the the half-domain |3^ , one 
corresponding to the Neumann case (N) and the other, 
with the additional minus sign, to the Dirichlet case (D): 



^00 (m) 



where A[f = A]^,'^"'^ 



(12) 
(13) 



ai—a2 = o, 



_ ,21 (m) 



ai— a2— a 



Furthermore the KKR determinant det Af°°(fc, a, r) 
dot Af {k, ai — a, a2 — a, r) of the full domain factorizes 
into the product of the determinants of these Neumann 
and Dirichlet KKR matrices 

00 

det M°°{k, a, r) ^ J| det Af°° (fc, a, r) 



ndet M°° (fc, a, r) det M°° (fc, a, r) 



det M°°{k, a, det Af°°(fc, a, r)!^ 



D 

(14) 



Thus the two-sphere system contains the Dirichlet 
sphere-plate system as special case, namely in the sym- 
metric limit a\ — a2 — a, the sphere-plate system is equal 
the Dirichlet case in the half-domain, and the pertinent 



Note that Figgis rotated by 90 degree relative to the conven- 
tions of Ref . |5Q | , such that our vertical symmetry plane is called 
"horizontal" there. 



multi-scattering determinant det M°l of the sphere-plate 
system is just given as 



det 



A/°l(fc,a,L°l) detAr°(fc,a,r^2(L°l+a)) 



(15) 

Note that the shortest surface-to-surface distance in the 
symmetric two-sphere case is given by L°° — r — 2a, 
whereas the shortest surface-to-surface distance in the 



sphere-plate case is i°l — hi 



see Fig.^ 



The exact expressions for the Casimir energy of the 
two-sphere Dirichlet problem, the symmetric two-sphere 
Dirichlet problem and the sphere-plate Dirichlet problem 
are given by the following integrals, respectively: 

£c'°^(ai,a2,i) 
he f°° 

= — / dfc4lndct Af°i°"(ifc4, 01,02, r = L+ai + a2) , 
27r Jo 

he r°° 

E^°{a,L) = — I dfc4lndetAf°°(ifc4,a,r = i + 2a), 
27r Jo 



f°l(a,L) 



he 
2^ 



poo 

I dfc4lndet Ai"°°(iA:4,a,r = 2(L-|-a))|j^ . 
Jo 



(16) 



In practice, these expressions have to be numerically inte- 
grated up to an upper value k™^^ which should be chosen 
large enough, such that the numerical value of the inte- 
gral is stable for some specified range of decimal places. 
For the sphere-plate case, k^'^^ ^ 10/L specifies a good 
choice. In order that the evaluation of the determinant 
of the matrix Mfi?{m]k,a^r = 2{L + a))\^ is stable, the 
upper value of k induces a maximal value Zniax for 1,1' 
and m (with \m\ <l,l' < ^max); namely 31, 46] 



''max _ 2 4 ' 



14a/L. 



(17) 



For small values of the separation L, the maximal angular 
momentum and therefore the size of the KKR matrices 
(which scale with P x I where the last factor results from 
the separable m quantum number) becomes rapidly very 
large and limits the range of applicability of this numer- 
ical computation of the exact integral to medium and 
large values of L, say, L > 0.1a, chiefly because of the 
handling of the 3j-symbols which scale with l'^. 



IV. LARGE-DISTANCE LIMIT 

As shown in Ref. (13|, it is possible to obtain signifi- 
cantly simpler expressions for the (integrated) density of 
states in the limit of very large separation or very small 
scatterers. If the wave length A = 27r/fc is much larger 
than the radii of the scatterers one can show that the 
KKR-matrix M{k) is given by 



[A/(fc)] 



33 



5" - {1-6" )Jj{k) 



exp{ikrjji) 



(18) 



' 33 



5 



vertical symmetry plane 



half-domain I 




half- 


-domain II 








a/ \ 




■« — >• 








J 









r 





FIG. 1: Two identical spheres of radius a at a center-to-center separation r. The vertical symmetry plane, the two half-domains 
and the surface-surface separation L°° in the global domain and 1/°' in one of the half-domains are shown. 



(see |51| for the analog in the 2D case) , where the indices we get the asymptotic result 

= 1, . . . ,N denote the scatterers, rjj' is the distance 
between their centers and fj{k) is the s-wave scattering jVf°°(™) 
amplitude on the j-th scatterer. 

In the case of two identical spheres of radius a, with 
their centers located at the distance r apart, one obtains 



{k,a,r) 



± 



hW(fca) 



V(2^i + l)(2/2 + l) exp(ifcr) ^ 

X Tl—^, OmO- (^O) 



ikr 



det M (k) = 1 cxp [ifc(2r - 2a)] 

= detM(A:)|jY detM(fc)|^ 



(19) 



^ife(r- 



-'^^l .(20) 



As usual, the determinant of this two-sphere system in 
the global domain factorizes into two sub-determinants 
calculated for one of the half-domains, one subject to 
Neumann boundary conditions and the other (with the 
minus sign) subject to Dirichlct boundary conditions. 

As mentioned, this expression can be derived from Eqs. 
(|12I13|) in the case of large center-to-center separations 
of spheres, fcr ^ 1: we can make use of the asymptotic 
expression for the spherical Hankel function 



(21) 



which is actually exact for li^^kr), such that the inverse 
multi-scattering matrix becomes 



ii,{ka) y/{2h + l)(2h + l) 



exp(ifcrO^ /Zi ^2 l\ fh h I 

^ ikr 2-.^^^ + ^) [m ~m O) [O 



(22) 



Using now the orthogonality relation for the 3j- 
symbols 



^(2/ + l) 



i=0 



mi m2 ml \m[ m!^ ^^J- Ormm[Om^m',, 



Since for fcr 3> 1 the only nontrivial matrix Mj°°J'''' is 
separable in li and I2, the corresponding determinant is 
simply given by 



det M°°{k,a,r) 

exp(2i/cr) 
(i/cr)^ 

exp(2ifcr) 



Ll=0 



hp'(fca) 



1 2 



ikr) 



,\2 



iA{ka)]^ 



(24) 



for the identical two-sphere case and by 

detM°l(fc,a,L) = det M°°{k,a,r = 2{L + a))\^ 
^^_exp^ - 



1=0 



hW(fca) 



^ 1 _ £^E(M (_,A(fca)) 



for the sphere-plate case. Here 



A(fta)^^(-l)'(2Z + l)-^ 



j;(fca) 



1=0 



hp^(fca) 



(25) 



(26) 



is the multipole expansion of the scattering amplitude. 
For ka ^ 1, the ratio ii{ka)/h[^\ka) becomes [s^ 

ji{ka) i(fca) 



h«(fca) 



12 X 32 X •• • X (2Z - 1)2 X {21 + 1) 

+ O {{kaf^+^+^'-'>) . (27) 



Thus the dominant effect comes from the s-wave scatter- 
ing only and the / > terms can be neglected. Conse- 
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quently, the scattering amplitude becomes 
. jo(fca) 



iA{ka) 



+ 0((fca^)) 



(28) 



= fcacxp(— ifca) 

which imphes that H24|l becomes H19|l and (|25|l becomes 
the Dirichlet part of (^01) . 

If the determinant in the global domain (|19() is inserted 
into the modified Krein formula © we get the following 
result for the integrated geometry-dependent part of the 
density of states in the s-wave limit of the two-sphere 
case [li (note r = L°° + 2a): 



s — wave V 



7r(L°°-|-2a)' 



■sm[2k{L°°+a)] + 0{{kaf 



(29) 

Analogously one gets the s-wave limit for the Dirich- 
let sphere-plate system by inserting the Dirichlet deter- 
minant of the half-domain, namely the second term of 
Eq.lEOI, into ®: 



s— wave 



k{2L°\ 



0{{ka 



(30) 



[note r = 2(L°I + a)]. 

Now, the Casimir energy for two identical Dirichlet 
spheres in the large L — L°° limit follows simply by in- 
serting A/'°°wavo: Eq. if^ . into the integral Q) and per- 
forming the Wick- rotation as in (jSJ: 



The large-distance scaling is therefore proportional to 
a/L^, in contrast to the Casimir-Polder energy between 
a molecule and a conducting plane which scales like 
a? /L"^ JH^. The difference is associated with the fact 
that the Casimir-Polder energy results from the induced 
dipole moment, whereas in the scalar scattering the 
monopole term gives the dominant contribution. In fact, 
if one omits the s-wave scattering term and starts in- 
stead with the the p-wave term in the scattering function 
A{ka), the scalar Casimir energy for the sphere-plate sys- 
tem would show a large-L behaviour 



wave 



Tra^ /90 



1 



(35) 



2L, 



which is compatible with the a?" jV^ scaling of the 
Casimir-Polder energy. Thus the correct large-distance 
behaviour of the scalar Casimir energy has nothin g to 
do with missing diffraction contributions (see Refs. [33, 
0,0,0) to the semiclassical trace formula as conjec- 
tured in Ref. and repeated in Ref. [2^. It is rather 
based on the replacement of the semiclassical summa- 
tion, which we will discuss in the next section and which 
is valid when many partial amplitudes contribute, by the 
leading term(s) in the multipole expansion |13ll5ll |. 

In summary, the asymptotic expressions for the 
Casimir energy are given by 



C700 

^s— wave 



(a,L) 



-hc- 



na 
'1? 



Ati{L- 
4 



a){L + 2aY 



where 



he TT 

167r2 90 



hcTT^ 

1440 

is the prefactor of the Casimir energy 



2 X 720 L3 



(31) 



(32) 



(33) 



of the corresponding scalar (Dirichlet) two-plate system 
where A is the area of the plates. Instead of performing 
the Wick rotation, one can compute these integrals along 
the real axis too. In this case one would have to include a 
convergence factor exp(— rye) and take the limit rj — )■ -1-0 
at the end of the calculation (in a similar manner to the 
Feynman prescription for propagators). 

Similarly, the Casimir energy for the Dirichlet sphere- 
plate case in the large L — L°l limit follows by inserting 

A/^i^avej Eq. (jSOJ, into the integral and performing 
the Wick-rotation: 



£ 



s — wave 



ia,L) 



-he 



L2 



47r(L + a)(2i + a) 
2 



1 



(34) 



£g°(L>a) - 4x 



90 KTi-a^ 



3.6958 



Kira 



(36) 



2L) 



in the symmetric two-sphere scalar Dirichlet case and by 



f°l(L » a) ^ 2 X -J « 1.8479-^ (37) 



in the sphere-plate scalar Dirichlet case. 



V. THE SEMICLASSICAL APPROXIMATION 

The semiclassical approximation of the scalar two- 
sphere problem in the framework of the Gutzwiller trace 
formula [2^ was pioneered in Refs. Here we will 

focus on the link between the scattering approach and 
these semiclassical methods. 

The sphere-plate system at surface-to-surface separa- 
tion L is a special case of the sphere-sphere case for two 
spheres of radii ai and a2 at center-to-center separation 
R = L + ai + a2 in the limit a2 — > oo. As shown in 
Ref. the integrated density of states for the two- 
sphere system follows semiclassically from the Gutzwiller 
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trace formula 126 



1 ^ exp(ni2fcL) 
= — im > 



n |det([M(ai,a2,L)]"-l)|i/2 



exp(ni2fcL) 



1 \ 
— Im y — 

TT |A+(ai, a2, i)" + A_(ai, 03, i)" - 2| 

E- 



sin(n2fcL) 



nTT |A+(ai, 02, L)" + A_(ai, aa, L)" - 2| 



, (38) 



where the periodic orbit is the bouncing-orbit between 
the spheres and the summation is over the repeats of 
this orbit. The matrix M(ai, a2, i) is the monodromy 
matrix of the sphere-sphere system and A+(ai,a2,i) = 
1/A_(ai, 02, L) is the double-degenerate leading eigen- 
value of this matrix, i.e.: 



A±(ai,a2,i) 



, 1 1 \ 
1 + 2L ( — + — +2 

ai a2 I aia2 



±\ 







\ 






+ 2 




02/ 


aia2/ 



-1, 

(39) 



which is identical with Eq. (3.11) of Ref. |17l |. 

In fact, the semiclassical expression (j38|l is consistent 
with the semiclassical limit to the exact expression (|^: 
In Rcf. 33] it was argued for the two-sphere case and in 
Ref. [31] it was shown for any A^-disk case that semiclas- 
sically 



det M{k) exp 



EE 

P 71—1 



p7r/2 



|det ([My 



1)1 



1/2 



where Mp and u 



(40) 

are the total geometrical length, 
the monodromy matrix and the Maslov index of the p-th 
primitive periodic orbit, respectively. The r.h.s. of H4U|) 
is the Gutzwiller-Voros zeta function [s^l ■ Note that for 
our scalar Dirichlet case there exists only one orbit, the 
bouncing-orbit for the two-sphere (two-disk) system with 
Ip — 2L, and that the Maslov index is simply I'p = 4 
because of the two Dirichlet reflections (per repeat). ^ 



^ For the asymptotic case kai ^ 1 the scattering amplitude 1261 of 
the previous section simplifies u nder the Debye approximation of 
the Bessel and Hankel function |52| and the replacement of the 
angular momentum sum by an integral; 

A{ka,) = f]i-l)'(2l + 1) if"'^ « i:^exp(-i2fca,). 

This implies 

2ifcr 

detM°l°2 ft! 1 ~ i-_A(fcai)A{fca2) » 1 - ^i2k(r-ai-a2) ^ 

(ikry Ar^ 



If the semiclassical expression (|38|l is inserted into the 
Casimir-energy integral (see (13j) 



f°° d 
" Jo dfc 

= -\hc / dfcAA,e(fc), 

Jo 



(41) 



one gets for the scalar sphere-sphere case, after a Wick 
rotation (as in the transition from Q to (|Hl): ® 

/•OO 

£°^°^{ai,a2,L) = -lhc / dfc AA°i°^ (ai, 03, A:) 
Jo 



00 



dfc4 exp(— n2A;4L) 



71 TT |A+(ai, 02, i)" + A_ (ai, 02, L)" - 2| 



he 



0102 



167r L^{ai+a2 + L) ^ 



00 ^ 
E 



1 



3 ai+a2 
he 0102 



3 aia2{ai+a2) 

4 

'1 



(n^-l) 



167r L2(ai + 02 + i) V90 
2L /5 1\ L{al + al) /5 1 



01+02 \7r 



0102(01+02) V"'^ 3 



(42) 



Here we applied the following identity 

A+(ai, 02, L)" + A_(ai, 02, i)" 
, _ _ L2 2 £2 

4n 





L 

H 


Ol 


02 






5 





- 2 

(n^ - 1) 



.(43) 



such that 



^ ^ asym 



Im In [1 



4r2 J 



sin [2A;(r — ai — 02)] 



The next term in the 1 / kr expansion of the Hankcl function I21i 
generates the correction 

« ^(l + ^ + ^)sin[2fe(r-ai-a2)] 



a\a2 



■ sin [2fc(r -- ai 



■ 0,2} 



47rr(r — a\ — 02) 

which is the n = 1 term of the Gutz wilier formula 1381 . consistent 
with the asymptotic limit L > ai ^ l/k. This is Eq. (13) of Il3l 
in the case ai =02. 
® The second line of Eq. 1421 , after the trivial integration over k^ , 
is identical with Eq. (3.20) of Ref. llTl if the latter is divided by a 
factor of two, as the zero modes are weighted there with a factor 
of 1 instead of 1/2. 
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which is exact for the case n = 1, and used 

CXD 4 oo 2 

El TT* , 1 

— — — and > — — — 
7^4 on /-^ n2 fi 



n=l 



90 



n=l 



6 



Particularly, in the case of two identical spheres of ra- 
dius a at center-to-center separation R = L-\-2a one gets 
a simplified expression 



he 



a 



167r L^{2a+L) 



L^i2a+L) 



1 



a 



a 

-'5 



(44) 



It is remarkable that the leading term is exactly equal 
to the plate-based prediction of the proximity-force ap- 
proximation for two identical spheres (ITJ (more details 
in Appendix El : 



coo 

Opiate PFA 



L^{2a + L) 



(45) 



As mentioned above, the sphere-plate system is a spe- 
cial case of the two-sphere system: 



he a 
na 



lim £°^°^{a,a2,L) 

a2 — *oo 
4 

' 'l , „ 

3 



1 



K 



i2 



1 



3/ a 



L 

- 4 
a 

0{[L/ar 



0{[L/ar) 



(46) 



This expression agrees with Eq. (11) of Rcf. 55] . Note 
that 



-S°\ia,L)<-K-. 



(47) 



Moreover, the leading rnra/L^ behaviour of H46() agrees 
with the leading terms of the plate-based and sphere- 
based proximity-force approximations for the Dirichlet 
sphere-plate problem, respectively (see Refs. [T9ll2^l2^ ) 



co\ _ 7ra 1 

''plate PFA ^ i 



L/a' 

'-■sphcrcPFA ~ ^2 ) 



(48) 



6^ 



l-(l + ^)ln(l + a/L) . (49) 



More details about the proximity-force approximation of 
the Dirichlet sphere-plate system can be found in Ap- 
pendix ^ 

The sphere-plate result (|46|) can also be derived from 
the Dirichlet part of the identical-two-sphere result using 



^ an 



Af°°\ 
-'^sc Id 



1 °° 1 



exp [in{r — 2a) fc] 



n |A+(a,r)" + A_(a,r)" - 2| 



(50) 



where (compare with Eqs. H38|l and (|39|l ') 




= lim A±(ai = a,a2,L = L°l) (51) 

a2 — ^oo 

under r = 2(L°I + a). Note that 

A±(a,r)2 = A±{ai = a,a2 ^ a,L ^ L°°) (52) 

for r = L°° + 2a, such that for two identical spheres the 
integrated density of states is semiclassically 

Jjoo ^ 1 jj^^ i exp (in2(r - 2a)k) 

t: ^^n|A+(a,r)2"-fA_(a,r)2»-2| ■ ^ ' 



VI. RESULTS 

The results and approximations, discussed in the pre- 
vious sections for the Casimir energy (a, L) for the 
scalar Dirichlet sphere-plate case in units of KTra/L^ = 
—hcTT^a/ {lAAQL'^) are shown in Fig.|21as function of the 
ratio L/a. This figure should be compared with Fig. 8 
of the world-line approach of Ref. and with Fig. 4 of 
the optical approach of Ref. f23j which both only present 
data for L/a < 4. The circles (a) represent the numeri- 
cally calculated exact expression H16f) for the sphere-plate 
system between L = 0.1a and L = 512a (the line is only 
there to guide the eye), the curve (b) shows the s-wave 
approximation H34|) . the line (c) represents the asymp- 
totic limit 1.847 of H37|l . the curve (d) represents the nu- 
merically calculated (Wick-rotated) integral (|41|l over the 
semiclassical expression (|50|l including all repeats, and 
the line (e) shows the analytical semiclassical formula 
(|46|l valid modulo 0(^{L/a)^) corrections. The curve (f) 
is the result of the plate-based proximity-force approx- 
imation H48|l . and the curve (g) represents the result of 
the sphere-based proximity- force approximation (|49|l . 

Our numerically calculated data (a) agree for L/a < 1 
with the in Ref. JJ,j published data of the world-line ap- 
proach, within the quoted (statistical) error bars which 
are already sizable at L/a = 1. We have included these 
world-line data, which cover the range from L = a/256 to 
L = 4a, in Fig[21and marked them by stars (h). Note that 
their central values are systematically on the low side in 
comparison to ours. All our data beyond L = 4a are pre- 
dictions. In the meantime, after the first version of this 
paper was submitted, there appeared in [s^ new data in 
the world-line approach with improved systematics which 
cover the extended range from L = a/512 to L = 16a and 
which have smaller, but still sizable error bars for L > a 
(see the crosses (k) in Fig.|2l. In the region of overlap, 
thus now also for the points L/a — 2,4,8,16, the new 
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world-line data do nicely agree with ours when the quoted 
statistical error bars are taken into account, although 
their central values are now systematically higher. It 
should also be remarked that our smallest point L = 0.1a 
is already affected by a sizable truncation error in the in- 
tegration. Of course this problem is a matter of numerics 
and not a matter of principle. 

Note that the s-wave approximation becomes a rea- 
sonable approximation to the exact data from L > 4a 
(it works very nicely for L > 15a in agreement with the 
estimate (|17|l ') and, moreover, that the exact data in- 
deed converge to the predicted asymptotic value ISO/tt"* 
of Eq. H37|l and do not show any Casimir-Polder a^/L'^ 
scaling. It is also interesting that the exact data, at least 
for the L values for which they could be calculated, are 
larger than 1 (in units of mra/L'^), whereas the semi- 
classical approximation and the proximity-force approx- 
imations are strictly less than 1 (in the same units). It 
should be noted that (at least) the upper error ranges 
of the old world-line data of Ref . > the complete new 
world-line data of Ref. [s^ (with the exception of the 
lower error ranges of the points L/a = 1/64, 1/128) and 
- below L ~ 0.1a - the results of the improved optical 
approach of Ref. are larger than 1 as well. 



The semiclassical calculation, starts out at a value 
1 (in units of Kna/L'^) for L/a —> 0, which was pre- 
dicted by the proximity- force a ppr oximation(s) |^ Q and 
confirmed in, e.g., Refs. [13 111 El IH E3- Then, 
for intermediate values of L/a, however, the semiclas- 
sical results, even though smaller than 1, are supe- 
rior to the results of the proximity-force approximation 
which become ambiguous 19, 22, 23]. For L ^ a, the 
contributions of the repeats of the bouncing-orbit are 
strongly suppressed and the numerically calculated semi- 
classical expression converges to the one-bounce result 
— ftca/(167rL^) = (90/7r*) x ktto/L^ which is smaller by 
a factor 1/2 than the exact asymptotic answer l(T7|l . 

VII. CONCLUSIONS 

We presented here an exact calculation of the scalar 
Casimir energy for the case of two spheres and a sphere 
and a plane. It is based on a new Krein-type formula 
which directly expresses the geometry-dependent part of 
the density of states by the inverse multi-scattering ma- 
trix of the pertinent scattering problem. Thus the corre- 
sponding Casimir energy follows from the energy-integral 
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over the multi-scattering phase shift (the logarithm of the 
multi-scattering determinant). The calculation is there- 
fore not plagued by subtractions of the single-sphere con- 
tributions or by a removal of diverging ultra-violet contri- 
butions. The asymptotic limit (|37|1 . the presented s-wave 
approximation (|34|l and all data with L > Aa arc totally 
new results. Moreover, contrary to claims in the liter- 
ature, the Casimir-Polder scaling of the scalar Casimir 
effect is excluded by our numerical and analytical calcu- 
lations. 

The two-sphere and sphere-plane cases are only two 
examples, and the formalism presented can be easily ex- 
tended to any number of spheres and planes as well (or 
disks and lines in two-dimensions). We have exemplified 
the calculation of the scalar Casimir energy only for the 
case of Dirichlet boundary conditions. One can replace 
the Dirichlet with Neumann boundary conditions or with 
any other conditions easily, or even replace the scatterers 
with arbitrary nonoverlapping potentials/nonideal reflec- 
tors. 

Aside from the exact results we have also discussed sev- 
eral approximation schemes, the large separation limit, 
the semiclassical limit and the proximity-force approxi- 
mation. The exact results (which are easy to calculate 
and are definitely simpler to evaluate than in a path inte- 
gral approach) should be looked upon as test examples for 
other approximate methods. These results already show 
that the proximity formula and the semiclassical/orbit 
approaches are limited to small separations only, typi- 
cally much smaller than the curvature radii of the two 
surfaces. 

One can make the argument that the dominant mo- 
menta of the fluctuating fields contributing to the 
Casimir energy at a separation L are of the order k 1/L 
(see Eq. (fT7|) ') and thus for large separations only the s- 
wave scattering is important. This is the main reason 
why the semiclassical approximation (which is valid when 
many partial amplitudes contribute) fails at large sepa- 
rations. 

Just the opposite is true at small separations, where 
semiclassics is pretty good and so is the proximity for- 
mula, when the separations are smaller by one or, re- 
spectively, two orders of magnitude than the curvature 
radii. The same type of analysis can be straightfor- 
wardly extended to cylinders, or even objects with less 
symmetry (in which case the corresponding individual 
T-matrices ap pearing in the inverse multi-scattering ma- 
trix (see Ref. |3lL 144. li^ ) become somewhat more com- 
plicated, due to the loss of spherical symmetry). 
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APPENDIX A: AMBIGUITY IN THE 
PROXIMITY-FORCE APPROXIMATION 



The proximity-force approximation (PFA) for the 
Casimir energy of two arbitrary smooth surfaces (with 
Dirichlet boundary conditions for the scalar-field case) is 
given by the surface integral over the Casimir energy per 
area which belongs to an equivalent parallel-plate system 
that locally follows the two surfaces H Sill IH HI] , 
i.e. 



SpFA ~ 



da e [z{a)] 



(Al) 



Here A is the area of one of the opposing surfaces which 
are locally separated by the (surface-dependent) distance 
z{a) and e [z{a)] is the corresponding Casimir energy per 
area. In general, the plate segment da is tangential to 
only one of the surfaces and therefore the local distance 
vector z(a) is perpendicular only to this surface and not 
to the other one. Thus Spfa is not uniquely defined, 
since the area A can be either one of the two opposing 
surfaces (or even one ficticious surface somewhere inbe- 
tween). Particularly, for the case of a sphere of radius 
a and a plate at shortest surface-to-surface separation 
L we get the following expression for the "sphere-based 
PFA" [lllllllll, where the local distance vector z(a) is 
perpendicular to the sphere. 



'"sphere PFA 



half —sphere 
7r/2 



|z(a)r 



-I 



Ka^27r / de sin{0) 



[L + a — acos(6')]'' 



r 2 I a \ a / 



l-(l + f)ln(l + a/L) 



(A2) 
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The coefficient k is again the prefactor k — — of the 
Casimir energy f "(i) = corresponding 
scalar (Dirichlet) two-plate system. 



On the other hand, the "plate-based PFA" |19ll22ll23| 
(i.e., the local distance vector z[a) is perpendicular to 
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the plate) follows from 



"--plate PFA 



K 



dx dy 



K I d(t> I dpp — 



z{x,y)\' 







na 







L + a — ^/a 



L2 1 + L/c 



(A3) 



Note that 



ira 



^sphere PFA < Opiate PFA < L'^ ' ^^^) 

Finally, the PFA for two spheres of a common radius 
a and shortest surface-to-surface separation L can be de- 
rived from the plate-based PFA ljA3|l of the sphere-plate 
case, with the fictitious plate on the vertical symmetry 
plane as in Fig.^ as follows: 



coo 

•Opiate PFA 



dx dy 3 

x2+a2<Q2 |2z(x,y)| 

2ir 



dpp- 



1 



ira 



8(|:)^ l + L/2a L^{2a + L) 



(A5) 



APPENDIX B: COMPARISON WITH THE 
TWO-DIMENSIONAL TWO-DISK AND 
DISK-LINE SYSTEMS 



where r = ri2 = t'21- The general two-disk system is 
characterized by a 6*21, symmetry. If the two disks have 
a common radius a, the global domain of this system 
is separated by the vertical symmetry axis (instead of 
plane) into two half-domains. The corresponding sym- 
metry group is now C2/1 and the KKR-matrix splits into 
two KKR matrices valid for the half-domains, one cor- 
responding to Neumann (N) boundary conditions of the 
scattering wave functions on the vertical symmetry axis, 
the other, with the additional minus sign, corresponding 
to the Dirichlet (D) case: 



M°° A 

-'"mm' In 

M°° I 

-'"mm' Id 



(B4) 
(B5) 



412 I 

^mm' |ai=a2=a 



I 

^mm' |ai=a2=a.- 



where A„„„/ 

The two-dimensional KKR-matrix in the large- 
distance limit reads 



[Mik)y 



(1 



<5«"')/f (e)^^^^^^ (B6) 



instead of (|18() . Here fj^{e) is the s-wave scattering 
amplitude in two dimensions. 

Since the asymptotic expression of the ordinary Hankel 
function reads 



" nkr 



(' 



TT TT 

1 I fcr — m 

2 4 



^ (-i)™H«(fcr), 
Eqs. (|B4IB5|) become asymptotically 



(B7) 



The two-dimensional ana log of the TV-sphere matrix in 
three-dimensions is [3lL l32| 



M: 



S'' S,n.-, + (1 - ) 



X (-l)'"'e'(™".'.-™'"-')H« ™,(fcr,,0, (Bl) 

where j,j' = 1,2, ... ,N are the labels of the N disks. 
The integers m, m! with —oo < m, m' < oo are the angu- 
lar momentum quantum numbers in two-dimensions, aj 
and usual, the radius of disk j and the dis- 

tance between the centers of disk j and j' , respectively. 
,h{kr) and }i[^\kr) are the ordinary Bessel and Hankel 
functions of first kind, and ajji is the angle of the of 
the ray from the origin of disk j to the one of disk j' as 
measured in the local coordinate system of disk j. 

The two-dimensional analog of the two-sphere KKR- 
type matrix ifinjl is pll4(il| 



^21 ^22 



_ ( Smm' A]^^, I ili2) 



^mm' "mm' 



with 



M°° (k n r) a_ / ^\mi ^mAka)_ 



X 1 



H«(fcr) (B8) 



instead of H23|l . This expression is separable in mi and 
m2. Therefore, the corresponding determinant is asymp- 
totically given by 



detM°°(fc,o,r) 



H«(fcr) 



E (-1) 



m— — oo 



, Jm(fca) 
il^m\ka) 



= 1- n^o^kr) [A{ka)]\ 

for the identical two-disk case and by 
det M°°{k,a,r) 



(B9) 



1 - h{,^^ (fcr) 



= 1 -H[,^^(fcr) [A{ka)] 



m— — oo 



. 3rnika) 
H^n^(fca) 



(BIO) 



Ail^ 



i-iy 



«/HW(fca 



) tt(i) n \ /-nn\ for the disk-line case. Note that the asymptotic rela- 
—^m-^m'ikr), (B3) ™ , ,ni., . , , . , 



tion H21|) is actually exact for hl^'^ (kr) and therefore holds 
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also for small fc- values. For the ordinary Hankel func- 
tion Hg^^(A;r), however, the corresponding formula l)B7|l 
is only an asymptotic relation. Since the Casimir en- 
ergy receives contributions from all k < 10/ L, it severely 

worsens the s-wave result if Hq^'' (fcr) were replaced by its 
asymptotic form from Hlj7|l . Here 



Jm(fca) 



(Bll) 



is the multipole expansion of the scattering amplitude in 
two dimensions. For fca ^ 1, the dominant effect comes 
from the m = (s-wave) contribution and the \m\ > 
terms can be neglected. The scattering amplitude be- 
comes Isil 



A{ka) = -^"i^"^ +0({ka)') 



0{{ka)'),{BU) 



where = 0.577 • • • is Euler's constant. For the asymp- 
totic case kr > kai ^ 1, one finds 



A{kai 



-2\ka 



m= — oo Hm^(fcai) 2 

for the scattering amplitude, such that 

2 g2ifer 



det M° 



1 



TT ifcr 



■A{kax)A{ka2) 



^ _ V^T^gi2fc(r-Qi-Q2) 



2r 



(B13) 



In fact, the two-dimensional analogs of the semiclassical 
expression H44|l for the identical two-sphere case and 146|l 
for the sphere-plate case are given by 



oo ^ 
5Z ri'i 



he 



^ n— 1 

167r^^^ L V L(l + L/2a) 



(B14) 



for the two-disk case and 
he 



E 



1 1 1 



2 7rn 2nL 2n\L 

n—l 



ftcC(3) 2 



= K2D 2- 



(B15) 



16^ L/L lVI 
for the disk-line case, where 

«2D = -T7^C(3) and C(3) = E " ^ 1-20205. (B16) 
Idtt ^-^ n'^ 

n—l 

The corresponding proximity-force approximation reads 
in the line-based scenario 



^linc PFA 



K2D / dx ■ 



1 



2 arctan , , 

2 \\/L(L+a) 



2a K2D- 



L+2a\ , y/LiL+2a) 



L{L + 2a)^jL{L + 2a) 



K2D — ~ K2D 2.22144^^. (B17) 



The proximity-force approximation for the disk-line sys- 
tem in the disk-based scenario is given by 



■^diskPFA - '^SDfl 



ir/2 



COS( 



«^2D 



-7r/2 (L -I- a — a COS ( 
K2D 2.22144 " 



\/2LVI 



(B18) 



as well. Note that in the limit L/a ^ the two- 
dimensional semiclassical approximation ljB15|) and the 
proximity-force approximations ljB17p and (jB18|) do ap- 
proximately agree, but are not identical. Furthermore, 
the exact result, in the range where it can be calculated, 
i.e., for L > 0.1a, does not scale as ^/aJTjjL, but rather 
as (ajLY^^ jL. The s-wave result is a good approxima- 
tion to the exact result for L > 10a, if }l'^\kr) is not 
replaced by its asymptotic form (|B7p . 
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